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Abstract —In this paper, we propose the Fourier frequency vec¬ 
tor (FFV), inherently, associated with multidimensional Fourier 
transform. With the help of FFV, we are able to provide physical 
meaning of so called negative frequencies in multidimensional 
Fourier transform (MDFT), which in turn provide multidi¬ 
mensional spatial and space-time series analysis. The complex 
exponential representation of sinusoidal function always yields 
two frequencies, negative frequency corresponding to positive fre¬ 
quency and vice versa, in the multidimensional Fourier spectrum. 
Thus, using the MDFT, we propose multidimensional Hilbert 
transform (MDHT) and associated multidimensional analytic 
signal (MDAS) with following properties: (a) the extra and 
redundant positive, negative, or both frequencies, introduced 
due to complex exponential representation of multidimensional 
Fourier spectrum, are suppressed, (b) real part of MDAS is 
original signal, (c) real and imaginary part of MDAS are 
orthogonal, and (d) the magnitude envelope of a original signal 
is obtained as the magnitude of its associated MDAS, which is 
the instantaneous amplitude of the MDAS. The proposed MDHT 
and associated DMAS are generalization of the ID HT and AS, 
respectively. We also provide the decomposition of an image into 
the AM-FM image model by the Fourier method and obtain 
explicit expression for the analytic image computation by 2D- 
DFT. 

Index Terms —The Fourier and Hilbert transform, analytic 
signal (AS), Fourier frequency vector (FFV), (time) frequency 
(rad/s) and spatial frequency (rad/m), single orthant Fourier 
transform (SOFT). 


I. Introduction 

The one-dimensional (ID) Hilbert Transform (HT) and 
analytic signal (AS) [l] have been used extensively in signal 
processing and information theory since their introduction. In 
this study, we use the Fourier theory to obtain the multidimen¬ 
sional HT (MDHT) and multidimensional AS (MDAS), as a 
generalization of ID cases, respectively. 

Let f(x) be a periodic function of time (or space), i.e. 
f(x) = f(x + T x ), \/x £ R, then the Fourier series (FS) 
of f(x) is defined as 


by 


a o 


dk 


bk 



f(x)dx, 

f(x) cos{ku>\x) dx, 
/(x) sin(fcu;ix) dx. 


k> 1 

k > 0. 


( 2 ) 


The Fourier amplitude spectrum for a real function is given 
by |To | = |ao|, |Tfc| = yjal + b\ and phase spectrum by 
ZFk = (f>k = tan ~ 1 (—bk/a,k), which implies that Fk can be 
written as 


To — ao, 


Fk = dk — jbk = 


2 

F, 



f{x)e { ~° k ^ x) dx, 

(3) 


i.e. Fk = |Tfc|e J0fc O d k = \F k \ cos(</> fc )^6 fc = |F fc | sin((?) fc ) 
and hence ([T]) can be written as /(x) = Y^kLo |Tfc| cos(fcwiX— 
$k )• 

In the Euler’s formula e W = cos (cf>) + j sin(0), the Hilbert 
showed that (with (j) = ut) the function sin(c6) is the HT 
of function cos(^). This yields the ±f degree phase-shift 
operator, which is a basic property of the HT. From the Euler’s 
formula, it is easy to show that cos ((f>) = \ [ e ^ + e - -^] and 
sin(</>) = F [e^ — e~^]. Using these results, equation (jTJ) can 
be written as 


OO 

f ( x ) = ex P(jkuJix) + c* k exp(-jkLOix)} <^> 

k—0 

oo 

f{x) = X ex PU ku} ix)] 

k=— oo 


(4) 


where Cq = do, c k = ( afc ^ bk ) for k ^ 0, from Q and (|4]i c k 
can be written as 


OO 

fix) = cos(/cwix) + bk sin(fcwix)] (1) 

k =0 

where loi = the coefficients do, a k and b k can be obtained 

A X 
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Ck = / f(x)exp(-jkuj 1 x)dx. (5) 

From Q, it is clear that the extra term, exp(— jkiO\x), which 
here corresponds to negative frequency, is introduced only 
due to complex exponential representation of sinusoids in <|TJ, 
otherwise this extra and redundant term is not required. 
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We consider the HT, for the sinusoidal functions, as the J 
degree phase delay operator, i.e. 

#{cosO)} = cos- ^) = sin(</>), 

n ( 6 ) 

H{sm[(p)} = sin(0 - -) = - cos(«/>), 

where <f> = cj^x for ID, = uj\X + for 2D and <j> = 
uj ix + - ■ ■ +ujmXm for MD function, the HT of function f(x) 
is denoted as H{f[x)} = fix). The HT of constant (ao) is 
zero, and it can be easily verified by 0 as H{ao cos(O)} = 
a o sin(0) = 0 and H{a 0 sin(7r/2)} = —a 0 cos(7r/2) = 0. 
The HT of complex sinusoidal (complex exponential) signal 
is multiplication of —j with signal itself (e.g. fT{exp(jwi)} = 
exp (j[uit — 7r/2]) = — j exp {j cot)). By using 0, we evaluate 
the Hilbert transform (HT) of (|T]» and obtain 

OO 

/ 0) = sin(fcwia;) — bk cos(fcwix)]. (7) 

k— o 

From 0 and 0 or simply from 0, the analytic signal (AS) 
is defined as @ 

OO 

z{x ) = 2 y^[c fc expjjk^x)], 

k =o 

z(x) = f(x)+jf(x) = r{x)e°^ x) 

where r(x) = [f 2 (x) + f 2 (x)]^ is a magnitude envelope of a 
signal f(x), 4>{x) = tan -1 [ real part of AS is original 
signal and imaginary part is the Hilbert transform of original 
signal. 

Observation: We observe that for real periodic function 

fix) 

OO 

y [a*; sin(Ax<;ix) — bj. cos(fco;ix)] = 0. (9) 

k——oo 

and this can be easily proved from 0 by equating imaginary 
part to zero. Notice the difference between 0 and 0, in the 
range of values of k in summation. 

Notice that this definition of AS satisfies the following 
properties: (a) the extra and redundant negative frequencies, 
introduced due to complex exponential representation, are 
suppressed i.e. k takes only positive values in 0, (b) real 
part of analytic signal is original signal, (c) real and imaginary 
part of AS are orthogonal, and (d) the magnitude envelope of 
a real signal is obtained as the magnitude of its associated AS, 
which is the instantaneous amplitude of the AS. 

The frequency, in Hertz (Hz), is defined as number of events 
or cycles per second. The fundamental period, T x , is the 
duration of one cycle, and is the reciprocal of the fundamental 
frequency , i.e. f\ = ^ ■ Hence, by definition, the frequency 
is positive physically. In real world, the negative frequency 
does not exist and the meaning of negative frequencies is only 
mathematical, and not physical. 

The extension of the HT and AS to the two-dimensional 
(2D) case and their applications to image processing has been 
limited, due to the non-uniqueness of the MDHT and MDAS. 
This fact has led to a variety of definitions with different 
approaches, to satisfy the ID conditions in the 2D case, such 


as the conventional 2D HT in the spatial domain 0 where 
the negative frequencies are not suppressed, directional HT 
and quaternionic 2D-AS |3), single orthant 2D-HT |[5J where 
real part of analytic signal is not original signal, monogenic 
signal that is based on the Riesz transform instead of the HT 
0, generalized radial HT (7), 2D-HT and its corresponding 
AS based on a combination of lD-HTs |8j. These approaches 
have some useful applications, however the 2D-AS obtained 
with them do not satisfy the all properties of ID AS. There 
are many interesting applications of the 2D-HT and the corre¬ 
sponding AS, such as edge detection |7], corner detection |4j, 
phase congruency calculations |9) and AM-FM image models 
G9’GD including others. 

The fundamental property of the analytic signal, especially 
from the viewpoint of image processing and recognition, is the 
split of identity |6). The complex AS in polar representation 
yields two local features, the instantaneous amplitude and 
instantaneous phase. These local features fulfill the property of 
invariance and equivariance ® i.e. the local phase depends 
only on the local structure, and the local amplitude depends 
only on the local energy (square of amplitude). If these 
local features are a complete description of a signal, they 
are said to perform a split of identity G3- The split of 
identity is valid only for band-limited signals with local zero 
mean property |)6]|. Thus, the split of identity is valid for 
all the sinusoidal functions and, hence, valid for the Fourier 
representation of a function. Therefore, the Fourier theory 
yields the AS representation of a signal that relies on an 
orthogonal decomposition of the structural information (local 
phase), and the energetic information (local amplitude). 

In this study, we present the multidimensional trigonometric 
Fourier series, generalize AS representation and HT using the 
Fourier theory for multidimensional function with following 
properties: (PI) In the Fourier spectrum of MDAS, the extra 
and redundant frequencies should be suppressed. (P2) Real 
part of MDAS is original signal. (P3) Real and imaginary part 
of MDAS are orthogonal, and (P4) the magnitude envelope of 
a MD real signal is obtained as the magnitude of its associated 
MDAS, which is the instantaneous amplitude of the MDAS. 

In the literature, there are various methods and applications 
0-1)201 of ID nonlinear and nonstationary time series. Re¬ 
cently, based on the Fourier theory, the Fourier decomposition 
method (FDM) for nonlinear and nonstationary time series 
analysis is proposed in |2J. We, in this study, present an 
extension of the FDM for (2D data) image signal, refer as 
2D-FDM, that yields multi-component AM-FM image model. 

This paper is organized as follows: the multidimensional 
Fourier series, Hilbert transform and analytic signal are dis¬ 
cussed in Section 2. The 2D-DTFT and associated analytic 
signal is discussed in Section 3. The 2D-FDM for AM-FM 
image model is discussed in Section 4. Section 5 introduces 
the single orthant MD-DTFT. Numerical results are given in 
Section 6. Section 7 presents conclusions. 

II. The multidimensional Fourier Series, Hilbert 

TRANSFORM AND ANALYTIC SIGNAL 

In this section, we discuss the 2D Fourier Series (2D-FS), 
MD-FS, introduce the concept of Fourier frequency vector 


SINGH P. el at.: ND-FT, ND-HT. 


3 


(FFV) and obtain associated analytic signal. Let /(x, y) be 
a periodic (with period T x ,T y ) function of 2D space (or ID 
space-time) i.e. f(x + T x ,y) = f(x,y + T y ) = f(x,y), 
\/x,y £ R, then the 2D Fourier series (2D-FS) of f(x,y ) 
can be defined as 


OO OO 

f(x,y) = ££[ cikj cos(kujiX + loj2v)+ 

1=0 k=0 

— 1 oo 

b k ,i sin(fcwix + lui 2 y) + ££[ a,k,i cos(fcwix + lu> 2 y) 


1= —oo k= 1 

+ b kj i sin(fcwix + lco 2 y) 


( 10 ) 


where u>i = ^ = 2 ttu, 1 x 2 = 5 ? = 2nv. Since, l is taking 
both negative and positive values, therefore individual spatial 
frequencies lu >2 can be positive as well as negative, but over all 
resultant spatial frequency (rad/m) is always positive, which is 
given by w = y/ ( kui \) 2 + ( lu > 2 ) 2 and negative sign only helps 
in determining the direction, 6 = tan _1 (fcwi/(tL> 2 ), of wave. 
Hence, with respect to Fourier theory, we refer them (fcuti, liv 2 ) 
as Fourier frequency vector (FFV) that has a magnitude as 
well as direction and can be written as u> = [few 1 IUI 2 ] (in 
physics, it is similar to a ‘wave vector’ in multidimensional 
systems). The coefficients 00 o> o k 1 and b k 1 can be obtained 
by 


0-0,0 — 


Ok,l — 


T x T y J_^J_^ 


T X Ty 


bk,i = 


T X Ty J_1 


f(x,y) dxd y, 

f(x, y) cos(fcwiX + luj 2 y) dx d y, 

f(x, y) sin(fcwix + luj 2 y) dx dy. 

( 11 ) 


The Fourier amplitude spec trum for a real function is given by 
l-Fo.ol = |a- 0,0 1^ \F k ,i\ = yja 2 k l + b 2 k l and phase spectrum by 
^F k ,i = 4> k ,i = tan _1 (—which implies that F k ,i 
can be written as 


J-y T x 

Fk,i = 7 =fr f 2 [ 2 f(x,y)e-^ x+l ^dxdy = 
-Lx-ly J-Im. 

o k ,i — jbk,i, F 0 ,o = ao,o- (12) 

We can also consider | [T()| with (a) limits of double sum from 
k = —00 to 00 , l = 0 to 00 and results would be same as 
there is no change in © (b) limits of double sum from 
k = —00 to 00 , l = —00 to 00 and values of a k 1 and b k ; in 
© would get divided by two. 

Discussion: There is a need for defining Fourier frequency 
vector in multidimensional Fourier representations. To demon¬ 
strate this need, let f{x,y) = cos(3wiX — du} 2 y) be a periodic 
signal (standing wave in 2D space, like image), with ui\ = 

!jj -2 7’ t = 3, T y = 4. From ( fTT| , we find only two 

solutions, a k ,i = 1 if (fe = 3, l = —4) or (k = —3, 
l = 4); bk,i = 0 for all fc, l. Interestingly, both solutions 
are same as cos(3wiX — 4uj 2 y) = cos(—3wiX + 4w2j/). In 


situation like this, we cannot avoid negative values of l (or 
fc), because there are no solutions with positive values (i.e. 
nonnegative integers No = {0,1,2, ••■}) of fc and l, which 
represent Fourier frequencies. This kind of situation does not 
arise in ID Fourier representation as fc takes values only 
zero onward in 0- Thus, for the multidimensional Fourier 
representation, we define the concept of FFV and in this case 
FFV is w = [3wi — 4 ^ 2 ] , which has a magnitude and 
direction, like any other vector. The magnitude of FFV is 
frequency w = w = yj (3wi) 2 + (4w 2 ) 2 , which is always 
positive by definition itself. The complex exponential represen¬ 
tation of this function, cos(3wiX — doj 2 y) = e J ( 3aJl x ~ 4 ^y) + 
e ~j (,3u 1 x-4,w 2 y), Produces the extra frequencies by second 
term e —j(3uix—4u2y) w 

Observation: We can also use ( [TO} for space-time (x,t) 
series analysis, e.g. ID wave equation, f(x,t) = cos(fcwif — 
lu> 2 x), where few 1 = = 2irkfi = 2nf, wave vector (or 

FFV) w = lu> 2 , wave number (or spatial frequency) |w| = 
^ and phase velocity v p = ^ = /A. 

By the Euler’s formula, we know that cos {<f>) = + 

e~^] and sin(</>) = — e~^], using these values we can 

write © as 


00 00 / , 1 \ 

/(*.») = ££[ ak ’ 1 k ' 1 exp[j(fcwix + hv 2 y)] 
1=0 k =0 

. ( a k,l +jb k i) ■ 

H- 2 -exp[-j(fcwix + lw 2 y) \ 


+ £ £ exp[ i(kuix + lu 2 y)\ 

I ——00 k =1 

. k,l jh,l) r -n 1 / M 
H--exp[—j (fcwix + ho 2 y)\ . 

This equation can be written as 


03) 


OO OO 

f(x,y) = ££[ Cfc.z exp[)(fcwix + lu 2 y)] 

1=0 k =0 

+ c l,i exp[-j(fcwix + lu} 2 y)\ 

— 1 OO 

+ £ £[ c k ,i exp[j'(fcwix + loj 2 y)\ 

1 =—00 k =1 


+ c k,i exp[-j(kujix + lw 2 y)] 




OO OO 

f(x,y) = E E Ck,i exp[j(fcwix + luj 2 y)\, 

k =—00 l =—00 


(14) 


where co,o = ao,o an d c k ,i = (akJ 2 jb ' i,,) for k,l = 

— 00 , • • • , 0,1, • • • , 00 . Hence, from © we can write 


T y T x 

Ck,i = 7r^r r / f(x,y)exp(-j[ku 1 x + lu} 2 y]) 

Lx J-y J-i± J-S IL 

dxdy. (15) 

From we observe that the extra term, exp[— j{ku>\X + 
lu> 2 y)\, which may corresponds to positive, negative, or both 
frequencies, is introduced only due to complex exponential 
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representation of sinusoids in ([TO]), otherwise this extra and 
redundant term is not required. 

By using (|5J, we evaluate the multidimensional Hilbert 
transform (MDHT) of ( [TO} and obtain 


OO OO 

/0, y ) = EE[ a fc ,z sin(fcwia; + loj 2 y) 

1=0 k —o 

— bk,i cos(ku>iX + lu> 2 y) 

— 1 OO 

+ E E[ ak,i sin(fccjia; + lui 2 y) 

l——oo k— 1 

— bk,i cos(ku)ix + lu 2 y) > 


(16) 


From © and ( fl6| i or simply from ( fl3| >, we define multidi¬ 
mensional AS (MDAS) 


z(x,y) = 2 EE Ck,i exp[j{ku)ix + lui 2 y)} 

1=0 k =0 
— 1 oo 

+ 2 EE Cfc,z exp[j(fcwia; + lu 2 y)\, 

l =—oo k=l 

z(X, y) = f(x, y ) + jf(x, y) = r(x, y)e 3 ^ x ’ v \ 


(17) 


where r(x,y) = [f 2 {x,y) + f 2 (x,y)] 2 is a magnitude enve¬ 
lope of a signal f(x,y), (j){x,y) = tanreal part 
of MDAS is original signal, imaginary part is the MDHT 
of original signal and MDAS has only originally present 
frequencies. Clearly, the complex exponential representation 
of sinusoidal function always yields two frequencies, negative 
frequency corresponding to positive frequency and vice versa, 
in the Fourier spectrum. Hence, MDAS suppress the extra and 
redundant positive, negative, or both frequencies, introduced 
due to complex exponential representation of multidimensional 
Fourier spectrum and satisfy all the MDAS properties PI to 
P4. Observation: We observe that for real periodic function 

f(x,y) 


E E [ ak,i sm(kujix+luj 2 y)—bk,i cos(kuiiX+luj 2 y)^ 

k =—00 l =—00 

= 0, (18) 

and this can be easily proved from ( [T4| by equating imaginary 
part to zero. 

The above discussion can be easily extended for MD-FS 
which can be written as 


OO OOOOOO 

f(xi,---,x M )= ■" E E E hi ,-am 

kM = —00 k^=—oo &2 —0 k\ —0 

cos(fciWia;i + • • • + 

+ b ku ... ,k M sia(kiu}ixi -)-h k M x> M XM)\ 

00 00 —1 00 

+ E ■" E EE 

kM =—00 ks=—oo /c2——00 k i = l 

cos(fciWia;i + • • • + k M x>MX M ) 

+ bk u - ,fc„ sin(fciwia:i H-h 1cmUmXm)\ > (19) 


fi x i) i x M ) y ' ■ ■ ■ y ' y ] c fei,-" ,kn 

/c M =-oo k 2 =—00 k 1=—00 

exp^fcituizi -I-h k M VMX M )], 


1 


c fci,••• ,fcjw — 


/(»!,••• i x m) 


T\ ■ ■ ■ T m J _ Zjl j 

x exp[—j(fc!WiXi + ■ ■ ■ + k M ^M x m )] dxt • • • dx M , (20) 


where c ku ... ,k M = (ak lt -,k M ~ jh u -,k M )/ 2 - We can easil y 
obtain MDHT by replacing cos with sin and sin with — cos 
in ( fl9| and obtain AS such that real part of AS is original 
signal, AS has positive resultant frequency u> = |u>| = 
\J (fciwi) 2 + • • • + (ijuUif) 2 and its FFV can be written as 
u> = [kiuii ■ ■ ■ fenfUAi] • The power of the FT can be 
realized from the fact that the HT and analytic representation 
of a signal are, inherently, present in the Fourier representation. 

Observation: We can also use © for (M — 1)D space- 
time series (xi,--- ,X(m- i),t) analysis, e.g. 2D wave equa¬ 
tion, f(x,y,t) = cos(kiojit — k 2 uj 2 x — k^ut^y), where 
fciwi = = 2nkifi = 2irf, wave vector (or FFV) 

l o = [k 2 u: 2 A: 3 W 3 ] T , wave number (or spatial frequency) 

|u>| = \J{k 2 oj 2 ) 2 + (fc 3 u; 3 ) 2 = ^ an d phase velocity v p = 

kiuix _ f\ 

M 

For non-periodic MD signal, f(x i,--- ,Xm), the MD 
Fourier transform (MDFT) and inverse MDFT (IMDFT) are 
defined as 


/ oo /»oo 

••• / fix 1,-" ,X M ) 

x exp[-j(wixi H-h %%)] dxi ■ ■ ■ dx M , 

/ oo /*oo 

■■■ 67[/i, • • • , /m] 

-OO J — OO 

exp[j(wixi -I-1- w m im)] d/i • • • d/M, (21) 

respectively. We define MDAS, for real signal f(x i, • • • , xm), 
as 


z(x i,---,x m ) = 2 ■■■ / C'f/i, - • • ,/m] 

— OO J — OO J 0 

exp[j(wixi -I-1- d/i • • • d/ M , (22) 

so its real part is original signal and imaginary part is HT of 
original signal. It is very convenient to take first frequency 
always positive (may be related to ID time or space) and rest 
of the frequencies (related to MD space) may be positive or 
negative or both, which is also clear from <[l), ( [T()| , ( fl9| ) and 
( |22| ), and this convention has been used throughout in this 
study. 

From the above discussions (ID equations G and G; 2 d 
equations ( flQ| and (fT7j)), we propose the following result: 

Proposition 2.1: Let s(x) be a real valued signal and 
z(x) = s(x) + j.s(x) is the AS representation of s(x). 
Then real-valued Fourier representation (RVFR) of s(x) and 
complex-valued Fourier representation (CVFR) of z(x) are 
same, i.e. S r (uj) = Z{lS), where S r {u>) is RVFR of s(x), 
Z(u) is CVFR of 2 (x), x = [xi • • • I'm] and u> = 
[oil • • • U) M \ ■ 
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For a 2D periodic signal, the RVFR is given by m and 
the CVFR is given by ( fl4] > and •dll- 

ill. The 2D-DTFT and associated analytic signal 

Let x[m, n] be a non-periodic and real function of time, then 
the 2D discrete time Fourier transform (2D-DTFT) of x[m, n] 
is defined as 

OO OO 

-Y(wi,w 2 ) = E E x[m, n] exp(— j[u>im + w 2 n]) 

m =—00 n =—00 

(23) 

and 2D inverse discrete time Fourier transform (2D-IDTFT) 
is debited as 

11 r r 

x[m,n} = — — J J X(wi,w 2 )exp(j[wiTO + w 2 n]) 

dwi dw 2 (24) 

It is easy to show that, from Eq. ( |23| , JY(— uii,— w 2 ) = 
X*(uj\ ,w 2 ), X(—cj i,cu 2 ) = X*(wi, —w 2 ). We rewrite 
Eq. ( |24l > as 


x[m, n\ = 

1 1 r f 


2 tt 27r . 


'0 Jo 


X(uii, w 2 ) exp(j[wim + w 2 n]) dwi dw 2 


r 0 


X(oji,ui2) exp(j[wim + w 2 n]) dcui dw 2 


' — TT J 0 
r 0 /*0 


+ 


X(uji, 0J2) exp(j[u;iTO + w 2 n]) dun du> 2 


' —7T «/ — 7T 
r 0 


-Y(un,u; 2 ) exp(j[wiTO + w 2 n]) dun dw 2 


' 0 j —IT 


(25) 


(26) 


In this Eq., hrst term (denoted as z \ [m, n] and 0 < un < 
7r,0 < w 2 < tt) is complex conjugate of third term (z*[m, n] 
and —7T < un < 0, —7r < w 2 < 0) and second term (denoted 
as z 2 [m, n] and —n < u > 2 < 0,0 < un < tt) is complex 
conjugate of fourth term (z| [m, n] and 0 < u > 2 < 7r, —tt < 
uj 1 < 0). As x[m,n] is real function, we can write 

x[m,n\ = Re{zi 4 [m,n]} = Re{z 2 o[m, n]} 

= Re{z 12 [m,n]} = f?e{z 34 [m, n]}, 

where f?e{.} denote real part of analytic signal (AS) 

Z14[m, n ] = 2 (zi[m, n] + z 2 [m, n]) = x[m , n] + jii 4 [m, n\, 
z 23 [m,n] = 2 (z[[m,n] + ^[m,™]) = x[m,n] + 3x23, [m, n ], 
0i 2 [m, n] = 2(zi[m, n] + z 2 [m, ?z]) = a?[m, n] + jxi 2 [m, n], 
z 34 [m, n] = 2(z*[m, n] + z 2 [m, ?t]) = :c[m, n] + j£ 34 [m, n], 

(27) 

where ii 4 [m,n] = —cc 23 [m.,n], ii 2 [m,n] = —i 24 [m,n] 

and subscripts denote the quadrants of the Fourier domain 
considered in AS representation. The DTFT of these AS, with 
their frequency supports, can be written as 

Z 14 (UJ 1 ,UJ 2 ), Wi <E [0,7r], U >2 G [—7T, 7r], 

Z 23 (W 1 ,W 2 ), WiG[-7 T,0], UJ 2 G [-7T,7r], 

Z 12 (uJl, UJ 2 ), UJ\ G [ — 7T,7r], W 2 G [0,7r], 
Z 34 (wi,w 2 ), Wi G [—7T, 7r], U> 2 G [-7T, 0]. 


Notice that with this dehnition of AS, the extra negative 
frequencies uji and u> 2 are suppressed in AS Z 44 [m, n] and 
z 12 [m, n], respectively, and this preserve the desired property 
for the AS. The AS s 23 [?n,n] and z 34 [m, n] are the complex 
conjugate of zi 4 [m,n] and zi 2 [to,7i], respectively. From the 
analytic image point of view, all four imaginary part of AS 
in yield different images. To obtain the the original 
signal we need to know the values of 2D-DTFT in hrst and 
second or third and fourth, hrst and fourth or second and third 
quadrants, which is also clear from the above discussions with 
equations ( |27| and ( |28[ i. It is to be noted that the only AS 
Zi 4 [to, n] coincides with the proposed AS in ( fT7) > and hence 
its imaginary part, /m{ 2 i 4 [m,n]}, coincides with proposed 
HT in (It)}. 

Example: Time-frequency analysis of 2D unit 
sample sequence. The unit sample sequence debned 
as S[m — mo,n — no] = 1 at m = mo,n = no 
and zero otherwise. We obtain the analytic 
representation of x[m,n] = S[m — mo,n — no] 4=> 
X{uji,uj 2 ) = exp(—j[wim 0 + w 2 n 0 ]) as 2 i 4 [m,n] 
= 2 2l [m,»l + 2«[m,„] = + 

* s s iv “ in 

its real part is original signal, i.e. x[m,n\ = 

[ •" i ;ri") )| ][ , "'irL;: l ,)) ] = h ™ -'»».»- h. 

obtain the phase of Z 44 as (f>[m,n\ = f(ni — mo) 
and hence wi[m, n] = which corresponds to half 
of the Nyquist frequency, and other frequency is 

w 2 [?n,n] = 0. Similarly, we obtain AS 2 i 2 [m,n] = 

2»k»] + 24[m,»] = + 

• rsin(7r(m——cos(7r(n—nn))) 1 • • /FT7T 

7 ——— 7 — ny x / —— [LLLL \, as given in d26D, 

J L 7 r(m—mo)ir{n—no) J’ © 1 i ’ 

its real part is original signal, i.e. x[m,n\ = 

1 [ •‘"ir~-^> >) i = J i m -”*»•»- "oi- w = 

obtain the phase of z 42 as <j>[m,n] = |(n — no) and hence 
w 2 [m,?t] = which corresponds to half of the Nyquist 
frequency, and other frequency is wi[m,n] = 0. 

If analytic signal z[m,n\ is debned as four times of brst 
quadrant, where 0 < w 4 < 7r,0 < w 2 < 7r, then we can 
observe that 

x[m,n\ Re{z[m,n]}, (29) 

where Re{.} denote real part of analytic signal (AS) 


r 

z\m. n\ = - 


/• TT /* TT 

/ / X(uj 1 ,uj2)exp(juj 1 m + juj 2 n) 

Jo Jo 

dwidaj 2 . (30) 

Notice that with this debnition of AS, the negative frequencies 
are suppressed but real part of AS is not original signal and 
this does not provide a phase delay of —tt/2 degree. 

IV. The 2D Fourier decomposition method and 
AM-FM image model 

We propose to use the following 2D Fourier decomposition 
method (2D-FDM) to obtain multi-component AM-FM image 
model debned as 

M 

f(x, y) = E gi ( x ' y ) + n ( x ’ y)> 

i= 1 


(28) 


( 31 ) 
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where n(x,y) is a noise representing any residue (constant or 
trend) components, and the gi(x,y) = cti(x, y) cos(cj>i(x, y)) 
are M monocomponent (zero mean narrow band) nonstationary 
image signals that represents the 2D Fourier intrinsic band 
functions (2D-FIBFs) for a image signal. We define 2D 
monocomponent signals that follow either of the following 
conditions 

<t>x(x,y) = > °> Vx > (32a) 

<t> v { x ,y) = d ^Qy V ^ > V ZF (32b) 

or both ( |32a| ) and ( |32b| ). 

Example: monocomponent 2D signal e.g. (a) g(x, y) = 
cos(3wia: — 4ui2 y) which follow ( |32a| ) or ( |32b| > (b) g(x, y) = 
cos(3wia: + 4w2j/) which follow both ( |32a| ) and ( |32b| i. 

We write analytic function n] from Eq. © as 


1 1 
7T 27T 


X(uii,U) 2 ) exp(j[tUiTO + W 2 tl]) dwi dui 2 = 


' — TV j 0 


M 

T eg [m, n] exp (jc/>i [ m , n \) (33) 

i =1 


where (with wio = 0, x\m = 7r) 


i i r r uii 

ai[m,n]exp(j<pi[m,n\) = -— / X(xi,x 2 ) 

exp(i[cuim + W2n]) dw! dtU2, (34) 

for i = 1, • • • , M. To obtain minimum number of AFIBFs in 
low to high frequency scan, for each i, start with uq( 7 ;_ij, 
increase and select the maximum value of u>u such that 
uti(j_i) < LJu < n and phase <pi[rn, n] is a monotonically 
increasing function with respect to (w.r.t.) u>u[m,n], i.e. 

a,i[m, n] > 0, uju[m, n) = ( cj>i[m + 1, n] - n]) >0, 

(35) 

for all m, n. Similarly, in high to low frequency scan, the 
lower and upper limits of integration in Eq. ( |34| will change 
to u>ii to wi(j_i), respectively, with wio = = 0, and 

we can obtain minimum number of AFIBFs by selecting the 
minimum value of tou such that 0 < uiu < and 

Eq. ( |35) > is satisfied. As we have represented Z\4 [m, n] by 
equations ([33]),([34j. l[35} and considered the decomposition 
w.r.t. only wi, the same way, we can represent ^i 2 [m,n] and 
consider decomposition w.r.t. only 102 - 

We proposed multivariate FDM (MFDM), based on zero- 
phase filtering (ZPF), in [|2j. Here, we propose to use MFDM 
for the decomposition of 2D signal (e.g. image) into AM-FM 
model defined in (J3T|. 


V. Single orthant MD-DTFT 

In this section, we consider multidimensional (2D and its 
extension to general MD) DTFT of real multidimensional 
signals. This approach evaluate DTFT of real signal in only 
first orthant and values in rest of the orthant are obtained by 
simple conjugation defined (for 2D and MD cases) as follows: 


A. Single orthant 2D-DTFT and 2D-analytic signal 

Let x[m, n) be a non-periodic and real function of time, then 
the 2D discrete time Fourier transform (2D-DTFT) of x[m, n ], 
we define as 

OO OO 

X(wi,cu 2 ) = E E x[m, n] exp(— iu>im — j^n) 

m =—00 n =—00 

(36) 

such that 

= i 2 = j 2 = -1, (37) 

where i and j are purely imaginary number and represent 
phase corresponding to frequencies oji and 022 , respectively. 
Since, there are two frequencies and W 2 and integration of 
frequency is phase, hence it is logical and natural to define 
two phases corresponding to two frequencies. Here, we define 
two conjugation, bar (—) and star (*) corresponding to i and 
j, respectively. Hence, we write (|36l> as 


X(ui, 0J2) — X r (u}i, W2) + ijXij(u>i,u>2) 

+ iXi(u>i, 012) + jXj(u>i, U2) ( 38 ) 

where, first term of the right side of this equation is first real 
part (FRP), second term is real part due to product of two 
imaginary parts corresponding to i and j, third and fourth 
term are imaginary parts. From ( [ 37 ] ) and (| 38 ]», we obtain the 
following equations for each of the four quadrants 

X\ (011,012) = [X r (u> 1,012) — Xij(u> 1,^2)] + 

j [Xi (wi ,u> 2) + Xj (wi ,0)2)], 

X2 (tut, CU2) = [X r {u>\, W2) + Xij(u>i, 012)] + 
j[—Xi(u>l,U)2) + Xj(cj 1,012)], 

X3 (uii , CU2) = [X r (u>i,U]2) — Xij (tUi, UI2)] — 
j [Xi (wi , 0J2) + Xj (oil, W2)], 

X4 (uii, CU2) = [X r (u> 4 ,U) 2 ) + Xjj(u>i,(jJ 2 )] — 
j[—Xi(w 1 , 0 ) 2 ) + Xj(u}i,U)2)]> 

where third quadrant is complex conjugate of first one, 
and fourth quadrant is complex conjugate of second one. It 
can be shown that the ( [ 23 ] ) and ( | 39 | ) are actually different 
representation of the same thing. We now explore the 
properties of this 2 D-DTFT: 

( 1 ) Symmetry: For real x[m,n], it is easy to show, 
from Eq. ([ 36 ], that (a) X(uti,uj2) = X(— wi, 1x2), (b) 
X* (wi, 1x2) = X{—lo\, —102), (c) X* (wi, UJ2) = X (oil, —W2). 
This clearly indicates that, if we know the value of X(u> 1,1x2) 
in first quadrant (0 < wi < n and 0 < 0J2 < tt) only, then we 
can obtain the value of X(101,1x2) in all the quadrants. The 
spectrum (|X(wi, W2)|) of this 2 D-DTFT is same as 2 D-DFT 
defined in ( | 23 | which is Hermitian symmetric. 

( 2 ) Periodicity: X(u>i,X2) is periodic in u>i, (02 with period 
27 t, i.e. X(x\,X2) = X(x\ + 2tt,X2 + 2 n). 

( 3 ) Shifting: x[m — mo,n — no] ^7 

exp(— im 0 uji) exp(—jn 0 X2)X(xi,x 2 )- 

(4) Modulation: exp(imwio) exp(jnw 2 o)a;[m, n] 

X(tOi — W10, X2 — X2o)- 

(5) Energy Conservation: ]C-oo S^°oo \x[m,n]\ 2 = 
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|X(wi,w 2 )| 2 dwidu; 2 , where |X(wi, w 2 )| 2 = 
X(wi,w 2 )X*(wi,u; 2 ). 

Now, it is strait forward to generalize the single orthant 
2D-DTFT case and obtain the single orthant MD-DTFT. Let 
x[ni,--- ,riM] be a non-periodic and real function of time, 
then the ND discrete time Fourier transform (ND-DTFT) of 
x[ni, ■ ■ ■ ,71m]- we define as 

oo oo 

X(wi, • • • ,w M ) = x[ni,---,n M \ 

n i=—oo um =—oo 

x exp(-jituini- Jm^mUm) (40) 

such that 

3k-Ji = -1, (41) 

where j-\ to jm are purely imaginary numbers and represent 
phase corresponding to frequencies uq to u>m, respectively. 
Since, there are M frequencies uq to ujm and integration of 
frequency is phase, hence it is logical and natural to define M 
phases corresponding to M frequencies. Here, we define M 
conjugations for each j\ to Jm- respectively. 

It is to be noted that all the above discussions have been 
done with MD-FS, MD-DTFT and associated MD-analytic 
signal. It is very easy to extend this discussion for all the 
variants of Hilbert transforms and Fourier theory including 
Fourier transform (FT), DTFS, DFT and FFT. 

VI. Numerical results and discussions 

To demonstrate that the HT and associated AS, defined 
through the FT in this paper, satisfy all the properties, we 
performed the following simulations and compare the results 
with some Hilbert transforms available in literature. 

A. Amplitude modulation, the Hilbert transform and Fourier 
representation 

Let f m (x,y) = A m cos(3uqa; — 4u> 2 y) be a modulating 
signal and f c (x,y) = A c cosfTOuqrr + 8u; 2 y) be a carrier 
signal. Figure FTI (middle), with uq = w 2 = §?-, T x = 3, 
T y = 4, A c = 1, A m = 1 and spatial sampling frequency 
u> s = 2nf s = 2n x 30 (rad/m). The AM modulated signal can 
be written as = [A m + fm(x, y)]A c cos(10uqa; + 

8 oj 2 y), Figure [T] (bottom). The envelope of this modulated 
signal is given by f env {x,y) = [A m + f m (x,y)}, which 
is shown in Figure [T| (top). Figure [2] shows the perfectly 
recovered AM signal (top), the HT of AM signal (middle) 
and envelope signal (bottom). 

B. Examples of 2D Analytic Signals 

Consider the two-dimensional harmonic signal f(xi,xf) = 
cos(uqaq) cos(w 2 a; 2 ). Using the Fourier series expansion 
of this signal is given by /(aq, x 2 ) = 5 cos(wia:i + w 2 at 2 ) + 
| cos(wi x-i — CO 2 X 2 ), which implies that f(x i,a: 2 ) signal is 
sum of two sinusoidal signals of phases (w \X\ + OJ 2 X 2 ) and 
(wi^i — 0 J 2 X 2 ), respectively. Using or ( p~6| ), and we 
obtain 2D-AS z(xi,x 2 ) = + \ e A“ 

This 2D-AS z{x\,xf) has amplitude dependent on x-[. x 2 and 



Fig. 1: The envelope signal f e nv(%, y ) = [1 V )] (top), carrier signal / c (#, y ) 

(middle) and AM signal f a m{ x , y) = fenv( x , y)fc( x , y) (bottom). 



Fig. 2: The recovered AM signal, f a m (x, y), from the FS (top), the HT of recovered AM 
signal /h(^ 5 y) (middle) and perfectly recovered envelope signal f e nv( x ,y ) (bottom). 


a nonlinear phase. This shows that the superposition of con¬ 
stant amplitude and linear phase signals result in a signal with 
variable amplitude and a nonlinear phase. Whereas, the 2D- 
AS with single-quadrant spectra |5| is given by V , ( a T> a; 2 ) = 
e i(u)ix 1 +u 2 x 2 ) shows f(x i,a; 2 ) has a constant amplitude 
and a linear phase in both coordinates. 

Other examples are presented in Table |T| and Table [II] for the 
comparative study of the proposed Hilbert transforms (HT) 
/(;i'i, .x' 2 ), which is based on the Fourier theory and phase 
delay, with the four Hilbert transforms |5l, namely the partial 
HT (PHT) in X\ direction f Xl (xi, x 2 ) with 2D-AS z(x 1, at 2 ) = 
f(xi,x 2 ) + jfx 1 (xi,x 2 ), PTH in x 2 direction f X2 (x i,x 2 ) 
with 2D-AS z(x i,x 2 ) = f(x i,x 2 ) + jf X2 {x i,* 2 ), total HT 
(THT) f T {xi,x 2 ) with 2D-AS z(xi,x 2 ) = f(x 1 ,x 2 ) + 
jfr{x 1,^2), and single-orthant HT (SOHT) fso(x i,x 2 ) with 
2D-AS z(x!,x 2 ) = f(x!,x 2 ) - Jt{xi,x 2 ) + j[f xi (x 1 ,x 2 ) + 
fx 2 ( x i,x 2 )]. From these examples it is clear that the only 
proposed HT with 2D-AS, presented in Table [II] is proving 
correct amplitude and phase in all the cases. 

Figurejijshows a signal /( x, y) and its proposed HT f(x, y), 
partial HT in x direction f x (x, y) and partial HT in y direction 
f y (x,y). This example clearly demonstrate that the proposed 
HT is not oriented to any direction, unlike PHT. 

From the Table [I] and Table [II] we observe that the proposed 
HT f{x,y) is closest to PTH in x direction f x (x,y). To 
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observe the difference between these two, we rewrite © as 

oo oo 

f(x,y) = 0-0,0 + ££[ ak,i cos(koj!X + luj 2 y) 

1=1 k= 1 

+ bk,i sin^uqrr + luj 2 y) + 

OO 

+ ^2 cos(kujix) + 6^,0 sin(fcu;ia;) 

k=1 
oo 

+ ^2 [a 0 ,z cos(luj 2 y) + bo,i sm(lLJ 2 y) 

i=i 

— 1 oo 

+ ££[ Ok,l COs(kuiiX + l(jJ 2 y) 


l = — OO /c—1 


+ b k ,i sin (ku)ix + luj 2 y) 


(42) 


TABLE III: An example signal x[m, n]. 


1 

-7 

30 

40 

5 

-6 

-70 

8 

-too 


TABLE IV: The DFT X[k, l] of an example signal x[m, n] of Table [15] 


-99 

6.00 - 71.01 i 

6 + 71.011 

85.50 - 174.07i 

54.00 + 27.7 li 

-91.50 - 139.43i 

85.50 + 174.07i 

-91.50 + 139.43i 

54.00 - 27.7li 


zero-phase filtering, to orthogonal 2D-FIBFs 1-7 in order 
of increasing frequency components by 2D-FDM. FIBF-1 is 
lowest frequency component and FIBF-7 is highest frequency 
component. Lena image is sum of FIBF-1 to FIBF-7. From 
Figure [TO} it is clear that this decomposition is not steered to 
any direction. 


and obtain a partial Hilbert transform in x direction as 

OO OO 

fx{x,y) =££ [ ak,i sin ( fewia; + luj 2 y) 

l -1 k —t 

- b k ,i cos(kuiX + luj 2 y) + 

OO 

+ ^2 Ko sin(/ca;i:r) — bk ,o cos(fcu;ia;) 

' L 

— 1 oo 

££[ a k ,i sin(/ctL>iir + lu 2 y) 


k —1 


+ 


E. Single orthant Fourier transform (SOFT) example 

To show how single orthant Fourier transform works, we 
take an example signal x[m,n\ which is shown in Table III 
The DFT of this signal is shown in Table |IV] The proposed 
single orthant DFT is shown in Table [V] where first orthant 
DFT (SODFT) elements are shown in bold numbers and values 
in rest of the orthant are obtained by simple conjugations. The 


l =—oo k=l 


elements of DFT Table IV can be obtained by the elements 
of SODFT Table [V| from the relation X[k,l] = (X r [k. /] + 
Xij[k,l ]) + j(Xi[k, l] + Xj[k,l]). 


— bk,i cos(ku)iX + lu 2 y) 


(43) 


Equations ( fT6| ) and ( |43| provide clear differences between the 
proposed HT and PHT in x direction. 


C. Edge detection 

The Hilbert transform can be seen as an edge detectors. 
The Figure [4] and [6] show examples of gray level checkerboard 
image and natural image, respectively. We have calculated the 
directional HT (DHT), corresponding to first and fourth quad¬ 
rant (DHT-1-4), first and second quadrant (DHT-1-2), using 
FFT based algorithm. It appears that these Hilbert transforms 
act as edge detection steered in the x and y directions. If we 
compare the two Hilbert components, we can see different kind 
of edges responding to the x and y directions. To observe the 
differences between the directional HT and proposed HT, we 
evaluated the proposed HT, AS, analytic phase and gradient 
of same gray level checkerboard image and natural image in 
Figure [5] and [7] respectively. 


D. 2D-FDM 

Lena image and its decomposition to orthogonal 2D FIBFs 
are shown in Figures [8] and [9] by considering first and fourth 
quadrant (w.r.t. wi), first and second quadrants (w.r.t. u> 2 ), 
respectively. It appears that these decompositions are steered 
in the x and y directions. If we compare these two decompo¬ 
sitions, we can see different kind of frequency components 


responding to the different, x and y, directions. Figure 10 


shows Lena image and its decomposition, using Fourier based 


VII. Conclusion 

In this study, we have proposed the Fourier frequency 
vector (FFV), inherently, associated with multidimensional 
Fourier transform. With the help of FFV, we have provided 
the physical meaning of so called negative frequencies in 
multidimensional Fourier transform (MDFT), which in turn 
provide multidimensional spatial and space-time series anal¬ 
ysis. We have shown that the complex exponential represen¬ 
tation of sinusoidal function always yields two frequencies, 
negative frequency corresponding to positive frequency and 
vice versa, in the multidimensional Fourier spectrum. Thus, 
using the MDFT, we have proposed multidimensional Hilbert 
transform (MDHT) and associated multidimensional analytic 
signal (MDAS) with following properties: (a) the extra and 
redundant positive, negative, or both frequencies, introduced 
due to complex exponential representation of multidimensional 
Fourier spectrum, are suppressed, (b) real part of MDAS is 
original multidimensional signal, (c) real and imaginary part 
of MDAS are orthogonal, and (d) the magnitude envelope of a 
original multidimensional signal is obtained as the magnitude 
of its associated MDAS, which is the instantaneous amplitude 
of the MDAS. We have also proposed the decomposition of an 
image into AM-FM image model by the Fourier method (2D- 
FDM) and obtain explicit expression for the analytic image 
computation by 2D-DFT. 
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TABLE I: The partial Hilbert transform (PHT) in x\ direction / X1 (#1 , #2), PTH in X 2 direction f X2 (#1 , X 2 ), total HT (THT) / t (# 1,2:2) and single-orthant HT (SOHT) 
fso(x 1,X 2 ). 


function f(x 1,^2) 

fx 1 (xi,X 2 ) 

fx 2 (x 1 ,X 2 ) 

f T (x 1,X 2 ) 

/so(*l.* 2 ) 

constant (ao) 

0 

0 

0 

0 

sin(cj 2 ^ 2 ) sin(cji^i) 

— sin(o; 2 ^ 2 ) cos(cjiiei) 

— cos(cd 2 ^ 2 ) sin(cjiaii) 

cos(cj 2 (r 2 ) cos(o;iaii) 

- sinjtdiXi + OJ2X2) 

cos(cj 2 ^ 2 ) sinfcjiici) 

— cos(a; 2 ai 2 ) cos(a;i:ri) 

sin(ci; 2 ^ 2 ) sin(cjifci) 

— sin(eJ 2 £ 2 ) cos(oaia;i) 

- cos(u;i 3 "i + 0120:2) 

sin(cD 2 ^ 2 ) cos(cjifci) 

sin(a; 2 ai 2 ) sin(a;i:ri) 

— cos (022^2) cos(cdiiei) 

— cos(cd 2 ^ 2 ) sin(cJirri) 

- cosjorixi + 0:20:2) 

cos(u;ia:i + UI2X2) 

sin(u;ia:i + UJ2X2) 

sinfotixi + 0:22:2) 

— cos(u;ixi +0J2X2) 

2 sin(o;ia:i + 0:20:2) 

cos(oJia;i — UJ2X2) 

sin(tJia:i — 0:2X2) 

— Sindian — UJ2X2) 

cos(ojia;i — 0:22:2) 

0 

sin(u;ia;i + U2X2) 

— cos(wiXi + U2X2) 

— cos(u;ia;i + 0:22:2) 

— sinfuriii +U2X2) 

— 2 cos(o:ixi +0:20:2) 

sin(ajia;i — U2X2) 

— cos(ojixi — OJ2X2) 

cos(u;ia;i — 0:22:2) 

sin(ujixi — 0:22:2) 

0 

COs{uJlXl) 

sin(cjifci) 

0 

0 

sinful#].) 

cos(oJ 20 ; 2 ) 

0 

sin(w 2 :E 2 ) 

0 

sin(o: 20 : 2 ) 

sin(o;irri) 

— cos(cuia;i) 

0 

0 

— cos(cjiaii) 

sin(uJ 2 £ 2 ) 

0 

— 003(0122:2) 

0 

— cos(oJ 2 X 2 ) 


TABLE II: The proposed Fourier theory and phase delay based 2D HT f(xx 2 ) and analytic signal (2D-AS) z(x 1 , . 7 : 2 )■ 


function f(x 1,2:2) 

proposed f(x i, #2) 

proposed 2 D-AS z(x 1,2:2) 

constant (ao) 

0 

ao 

sin(cD 2 ^ 2 ) sin(cjiaii) 

— sin(cj 2 ^ 2 ) cos(coi:ri) 

i[ e j(o:i3:1-0:22:2) _ e 0 '(o> 1 xi+u: 2 a: 2 )] 

cos(o; 2 ^ 2 ) sin(cjia;i) 

— cos(co 2 ^ 2 ) cos(cjiaii) 

i [gjOuiXi- 0 ) 2 * 2 -f ) _|_ e 0 '(o> 1 xi+u: 2 a: 2 -f )l 

sin(o; 2 ^ 2 ) cos(cjiaii) 

sin(o; 2 ^ 2 ) sin(o;ia:i) 

I [+>(01x1-02x2 + 5 ) _|_ £ j (01 xx +o 2 x 2 - f ) j 

cos(o:io:i + 0:20:2) 

sin(oria:i + 0:20:2) 

e j(ui 1 x 1 +ui 2 x 2 ) 

cos(o:io;i — 0:20:2) 

sin(o:io:i — 0:20:2) 

e j(ui 1 x 1 -ix 2 x 2 ) 

sin(o:io;i + o: 2 x 2 ) 

— cos(oiio:i + 0:20:2) 

e j(u>ixi+ui 2 x 2 -^) 

sin(o:io:i — 0:20:2) 

— cos(oiio:i - 0:20:2) 

e j(tx 1 x 1 -ui 2 x 2 -^) 

cos(o;iiei) 

sin(cjiaii) 

gj(01XI) 

003(0:22:2) 

sin(o: 20 : 2 ) 

e 0(o2x2) 

sin(o;ia:i) 

— cos(o;iiri) 

e y(oia:i-f ) 

sin(o; 2 o: 2 ) 

— cos(o> 2 o:2) 

gf(" 2 X 2 -f ) 


Original signal f(x,y) 



Proposed HT f(x, y) 




Partial HT in y direction j y (x. y) 



Fig. 3: Original signal /(a;, y ) and its proposed HT f(x, y). partial HT in x direction (a:, y) and partial HT in y direction f y fa;, y ). 


TABLE V: The proposed single orthant DFT elements, of example signal a:[m, n] of Table |lll| in first orthant are shown in bold numbers and values in rest of the orthant are 
obtained by simple conjugations. 


A r 1 A i j , A i , J(. j 

-A r j A ij , -X i , JC j 

Xr.Xi^X^Xj 

- 99 , 0 , 0 , 0 

6 . 0000 , 0 , 0 , - 71.0141 

6.0000, 0, 0, 71.0141 

85 . 5000 , 0 , - 174 . 0711 , 0 

- 18 . 7500 , 72 . 7500 , - 55 . 8586 , 83.5715 

-18.7500, -72.7500, -55.8586, -83.5715 

85.5000, 0, 174.0711, 0 

-18.7500, -72.7500, 55.8586, 83.5715 

-18.7500, 72.7500, 55.8586, -83.5715 


Appendix 

The 2D-DFT and analytic image 
In this appendix, we derive analytic image by 2D-DFT. 
Let x[m,n\ be a real function, then the 2D discrete Fourier 


transform (2D-DFT) of x[m,n] is defined as 

M-l N~1 


X(k, l) = 


1 


MN 


EE 

m= 0n—0 


x\m, n\e 


(44) 
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Fig. 4: Directional HT (DHT) (a) Checkerboard image x[m,n] (b) DHT x \4 [m,n\ (c) DHT X 23 [m,n] (d) DHT &12 [m,n\ (e) DHT £34 [m,n\ (f) analytic amplitude 
\z\&[m, n]| (g) analytic amplitude |zi 2 [m, n]| (h) analytic phase Z.zi 4 [m , n] (AP-1-4) (i) analytic phase Z.zi 2 [m, n] (AP-1-2) (j) gradient 0m .[m, n] of AP-1-4 (k) gradient 
4> n [m,n\ of AP-1-4 (1) \f 0^ [m, n] + 0^[m,n] of quadrant 1-4 (m) gradient 0 m [m,n] of AP-1-2 (n) gradient 0 ri [m,n] of AP-1-2 (o) \J 0^ [m, n] + 0^ [m, n] of 
quadrant 1 - 2 . 



Fig. 5: Proposed HT (a) Checkerboard image x[m,n\ (b) HT £± 4 [m, n] (c) HT £ 23 [m,n\ (d) analytic amplitude |.zi 4 [m,n]| (e) analytic phase Z.zi 4 [m, n\ (AP-1-4) (f) 
gradient 0 m [m, n] of AP-1-4 (g) gradient 0 ri [m, n] of AP-1-4 (h) \J 0^ [m, n] + 0^ [m, n] of quadrant 1-4. 


and 2D inverse discrete Fourier transform (2D-IDFT) is de¬ 
fined as 

M—l N—l 

x[m,n] = ^ ^ X[k, (45) 

fc—o (=0 

In ( |44| . for odd numbers (M = N ), there is only 
one real term X[0,0] and (MN — l)/2 terms are com¬ 
plex conjugate of the rest (AIN — l)/2 terms; and for 
even numbers (M = N), there are only four real terms 


X[0,0],X[0,N/2],X[M/2,0],X[M/2,N/2] and (MN - 
4)/2 terms are complex conjugate of the rest (MN — 4) /2 
terms. From the above discussions and using the conjugate 
symmetry of 2D-DFT, we obtain 2D-AS for odd numbers 


















Fig. 6 : Directional HT (DHT) (a) Lena image x[m,n\ (b) DHT #14 [m, n] (c) DHT #23 [m,n\ (d) DHT £12 [tti, 71 ] (e) DHT £ 34 ( 771 , 71 ] (f) analytic amplitude 1 214 ( 771 , 71 ] | 
(g) analytic amplitude 1 212 ( 771 , n]| (h) analytic phase Zzi4[77i, n] (AP-1-4) (i) analytic phase Z.zi 2 [m , 71 ] (AP-1-2) (j) gradient 4> rn [m, 71 ] of AP-1-4 (k) gradient (p n [m, n] of 
AP- 1-4 (1) [ 771 , 71 ] + 0^[t 71, 7i] of quadrant 1-4 (m) gradient 0 m [r7i, 71 ] of AP-1-2 (n) gradient <p n [m, n] of AP-1-2 (o) yj [tti, n] + [tti, n] of quadrant 1-2. 



Fig. 7: Proposed HT (a) Checkerboard image #[771,71] (b) HT #14(771,71] (c) HT #23(771,71] (d) analytic amplitude I214[tti, ti] | (e) analytic phase Z214 [m, n] (AP-1-4) (f) 
gradient 0 m. [tti, 71] of AP-1-4 (g) gradient 0 ^ [tti, ti] of AP-1-4 (h) yj 0 ^ [tti, ti] + [tti, ti] of quadrant 1-4. 


(M = N) as 


(M—1)/2 (AT-l)/2 

zu[m, n] = 2 X[k, l]e^^ +w ^ 

k —0 1=0 

(M—1) (JV—1)/2 

+ 2 ^ -X(0,0) (46) 

A—(A/ • J)/2 i=l 


and for even numbers (M = N) as 

M/2 N/2 

Zi4[m, n] = 2 ^ ^ X[fc, 
fc=0 A=0 
M-l N/2— 1 

+ 2 ^ ^ X[k,l]e i ^ +l ^- 

k=M/ 2+1 Z=1 

X[0, 0] - X[0, N/2} - X[M/ 2,0] - X[M/2, iV/2], (47) 
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Lena image 
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FIBF -3 
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FIBF-6 
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Fig. 8: Lena image and its decomposition, w.r.t. cj i by p4| , to orthogonal 2D-FIBFs 1-7 in order of increasing frequency components by 2D-FDM. FIBF-1 is lowest component 
and FIBF-7 is highest frequency component. Lena image is sum of all FIBFs 1-7. 



Fig. 9: Lena image and its decomposition, w.r.t. 0 J 2 , to orthogonal 2D-FIBFs 1-7 in increasing frequency components by 2D-FDM. FIBF-1 is lowest frequency component and 
FIBF-7 is highest frequency component. Lena image is sum of all FIBFs 1-7. 


where real part of AS is original signal, i.e 
x[m, n] = Re{zu[m , n ]} and imaginary part of AS is the HT 
of original signal, i.e. x[m,n] = Im{zu[m, n]}. This 2D-AS 
has been obtained by considering the first and fourth quadrants 
of 2D DFT, the second and third quadrants of 2D DFT. 
These 2D-AS (analytic image) computation and 2D Hilbert 
transform can be easily implemented with 2D-FFT algorithms, 
e.g. MATLAB implementation is presented in Algorithm A. 


Algorithm A: MATLAB code for the proposed analytic 
image and 2D Hilbert transform computation by 2D-FFT. 
Case 1: when M and N are even numbers. 
X=fft2(x);[M,N]=size(x) % compute 2D-FFT and size. 

% create mask. 

Xm=zeros(M,N); Xm(l:((M/2)+l),l:((N/2)+l))=2; 
Xm(((M/2)+2):M,2:(N/2))=2;% take 1st & 4th quadrant. 
Xm(l,l)=l; Xm(l,N/2+l)=l; Xm(M/2+l,l)=l; 

Xm(M/2+1 ,N/2+1 )=1; 

% mask 2D-FFT and compute inverse to obtain AS. 
tmp=Xm. *X; zl4=ifft2(tmp); 






























SINGH P. el at.: ND-FT, ND-HT. 
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Lena image 




FIBF-4 



FIBF-5 



FIBF-6 



FIBF-7 



Fig. 10: Lena image and its decomposition, using Fourier based zero-phase filtering (ZPF), to orthogonal 2D-FIBFs 1-7 in order of increasing frequency components by 2D-FDM. 
FIBF-1 is lowest frequency component and FIBF-7 is highest frequency component. Lena image is sum of all FIBFs 1-7. 


Case 2: when M and N are odd numbers. 

X=fft2(x); [M,NJ=size(x) % compute 2D-FFT and size. 

% create mask. 

Xm=zeros(M,N); Xm(l:(M+l)/2,l:(N+l)/2)=2; Xm(l,l)=l; 
Xm((M+3)/2:M,2:(N+l)/2)=2; % take 1st & 4th quadrant. 
% mask 2D-FFT and compute inverse to obtain AS. 
tmp=Xm. *X; zl4=ifft2(tmp); 
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